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A systematic study of the /n'^ ratio in heavy-ion collisions with the same neutron/proton 
ratio but different masses can help single out effects of the nuclear mean field on pion production. 
Based on simulations using the IBUU04 transport model, it is found that the 7r^/7r+ ratio in 
head-on collisions of ****Ca-|-''**Ca, ^^''Sn -l-^^^'Sn and ^''^Au+^^'^Au at beam energies from 0.25 to 
0.6 GeV/nucleon increases with increasing the system size or decreasing the beam energies. A 
comprehensive analysis of the dynamical isospin fractionation and the tt' /n^ ratio as well as their 
time evolution and spatial distributions demonstrates clearly that the -k^ /-k^ ratio is an effective 
probe of the high-density behavior of the nuclear symmetry energy. 

PACS numbers: 25.70.-z, 25.60.-t, 25.80.Ls, 24.10. Lx 



The high-density (HD) behavior of the nuclear symme- 
try energy Esymip) has long been regarded as among the 
most uncertain properties of dense neutron-rich nuclear 
matter P, 0, H, 01 ■ It is very essential for understand- 
ing not only many fundamental astrophysical phenom- 
ena 0, 0, 0] but also novel features of high energy heavy 
ion reactions especially those induced by rare isotopes 
0, i, i, 0, [Hi- The 7r-/7r+ ratio in heavy-ion colli- 
sions has been known as a particularly sensitive probe 
of the HD behavior of the Esy^{p) Based on the 
IBUU04 transport model analysis of the 7r~/7r+ data 
from the FOPI collaboration [IjI, circumstantial evi- 
dence suggesting a rather soft -E'sym(p) compared to the 
widely used APR ( Akmal-Pandharipande-Ravenhall) 
prediction[3l was reported very recently [l^. The 
IBUU04 transport model used the single particle poten- 
tial U corresponding to a modified Gogny Momentum 
Dependent Interaction (MDI) [l^, i.e., for a nucleon 
with momentum p and isospin r. 
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In the above equation the isospin r = 1/2 (—1/2) for 
neutrons (protons). The coefficients Au{x) and Ai{x) 
are [itI 



Au{x) 



-95.98 -X- 



2B 



Ai{x) 



-120.57 -f-x- 



2B 



1-1 
(2) 



'Electronic address: | xiaozg@tsinghua.edu.cn] 



4/3, B = 106.35 
103.40 MeV and 



The values of the parameters are a 
MeV, Crr = -11.70 MeV, Cr -r = 

A " ' 

at the saturation density po UM ■ ^'^^ asymmetric nuclear 
matter at zero temperature, the MDI symmetry energy 
can be written as 



which is the Fermi momentum of nuclear matter 
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where Pf = ft(37r2|)^/3 is the Fermi momentum for sym- 
metric nuclear matter at density p. We note here that 
since the Ai{x) - Au{x) = -24.59 -f ABx/{cr -I- 1) ac- 
cording to Eq. ([2]), the Esym{p) depends linearly on the 
parameter a; at a given density except at po where the 
i?sym(po) is fixed at 30.54 MeV. Shown in the inset of 
Fig. 1 are examples of the Esym{p) with x=l. and -2, 
respectively. The strengthes of the corresponding isovec- 
tor (symmetry) potential estimated from ([/„ — Up)/25 
are shown for 3 typical densities in the main frame of 
Fig. [T] At the normal density po, by design, the symme- 
try potential is independent of x and is consistent with 
the Lane potential extracted from the experimental data 
of nucleon-nucleus scatterings and the (p,n) charge ex- 
change reactions It is necessary to emphasize that 
it is the symmetry potential, not the symmetry energy, 
that enters as a direct input in all transport model sim- 
ulations. It is thus important to point out the key char- 
acteristics of the symmetry potentials. With x=0, the 
symmetry potential is weak but remains mostly posi- 
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tive at all densities in the momentum range considered. 
With x=-2 leading to the stifFer symmetry energy, the 
symmetry potential is positive at all densities. With x=l 
leading to the softer symmetry energy, however, the sym- 
metry potential is negative at supra-saturation densities. 

It was found that only the very soft -Bsym(p) with 
X = 1 can well reproduce the FOPI data . The stiffer 
Esym (p) with X = (which resembles very well the APR 
prediction up to about 3.5po [ll]) under-predicts signifi- 
cantly the data as in earlier IQMD [l^ calculations 
using an i?sym(p) very similar to the APR prediction or 
the MDI Esymip) with x = 0. The interesting ramifi- 
cations in both nuclear physics and astrophysics of this 
finding about the HD behavior of the i?sym(p) strongly 
calls for additional theoretical studies and experimental 
tests. Since pions always undergo strong final state in- 
teractions, the latter may distort the information carried 
by the ■k~ j-K^ messenger about the HD i?sym(p) • In this 
work, based on the IBUU04 transport model simulations 
of heavy-ion collisions with the same neutron/proton ra- 
tio but different masses we further investigate how reli- 
able the 7r~/7r+ ratio is in probing the HD behavior of 
the Esymip) ■ 



> 



to 



700 
600 
500 
400 
300 




200 - 



i3 100 







-100 



-200 



P/Po 



k(fm"^) 



FIG. 1: (Color online) The symmetry potential as a function 
of nucleon momentum. With p = po, the curves are identical 
for x=-2, and 1. The inset shows the density dependence 
of the symmetry energy with x—-2, and 1, respectively. 

To understand the advantage of comparing system- 
atically the tt^/tt"*" ratios in heavy- ion reactions with 
the same neutron/proton ratio but different masses, we 
start by recalling the first-chance nucleon-nucleon col- 
lision model for pion production and the dynamical 
isospin fractionation mechanism for enhancing the neu- 
tron/proton ratio of the HD region and thus the tt~ /tt'^ 



ratio from there. Near the pion production threshold 
of about 300 MeV, it is reasonable to assume that only 
first chance inelastic nucleon-nucleon collisions produce 
pions. Since the 7r~ and tt"*" are mostly produced from 
nn np + 7r~ and pp np + tt"*", respectively, and 
the pn collisions contribute equally to the production 
of both TT^ and 7r+. the primordial 7r~/7r"*" ratio is ex- 
pected to be proportional to the neutron/proton ratio 
squared, i.e., [N / Z)'^^^^^, of the participant region. In- 
deed, more detailed estimate based on the A(1232) iso- 
bar model [lOl predicts a primordial tt^/tt^ ratio of 
i?i,ob = (7r-/^+),cs = (5iV' + NZ)/{5Z^ + NZ) « 
{N/Z)\^^^^. Due to the detailed balance, the reabsorp- 
tion of 7r~ and tt"*" mainly through the A resonances 
correlate proportionally to their production, they are 
thus not expected to change significantly the primor- 
dial 7T^/7r+ ratio in the absence of possible in-medium 
effects of the A resonances (see, however, refs. [1^ and 
(soj ) . The Coulomb force in the final state is expected to 
affect the differential 7r^/7r+ ratio as a function of pion 
momentum but not the integrated 7r~/7r+ ratio. The 
final 7r~/7r+ multiplicity ratio is thus a direct measure 
of the isospin asymmetry {N/Z)({cnsc in the early stage 
of the reaction. While the latter is determined by the 
HD Esymip) through the dynamical isospin fractiona- 
tion phenomenon [2l|, [13, HI, [13, [1E| , namely, the par- 
ticipant region is more neutron-rich (poor) if the value 
of the HD Esymip) there is lower (higher). In terms of 
the reaction dynamics for pion production within trans- 
port models, the Tr~ /tt~^ ratio depends on the (Af/Z)ciense 
which is determined by the isovector part of the nuclear 
mean-field. Effects of the latter depends on the den- 
sity gradients reached and the duration of the reaction. 
The density range and reaction time can be varied by 
varying the masses of the reaction system and the beam 
energy. In doing so, to examine the effects of the HD 
Esymip) , if is better to use reactions of the same neu- 
tron/proton ratio so that the initial isospin asymmetry 
of the reaction system remains the same and thus allows 
us to examine clearly effects of the isospin fractionation. 
The observed effects can then be essentially attributed 
to the variation of the dynamical isospin fractionation 
due to the changing density and reaction time. In the 
following, we compare head-on reactions of *^Ca-|-^*Ca, 
^^^Sn-|-^^^Sn and ^^^Au4-"'^^^Au at beam energies from 
0.25 to 0.6 GeV/nucleon . Their isospin asymmetries 
are approximately the same, more quantitatively, 1.40, 
1.48 and 1.49, respectively. 

Shown in Figure [2] are the excitation functions of the 
7r~/7r+ ratio with the £'sym(p) of a; = 1 and x = 0, re- 
spectively. The error bars presented in the plots are 
of only statistical origin. The 7r~/7r+ ratio increases 
with decreasing the beam energy and exceeds the iso- 
bar model prediction indicated by the dashed line in 
the plots. This trend is again consistent with the FOPI 
data in Au-|-Au reactions and other model calculations 
[Tsj . In the energy regime considered here, pions are 
mainly produced through the A isobar. If we denote 
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FIG. 2: (Color online) Excitation function of the 7r~/7r^ ratio 
using the MDI interaction with x = 1 (open square) and 
(open triangle) for the head-on collisions of *^Ca-|-**Ca, ^'^''Sn 
+i2*Sn and "^Au+^^^Au, respectively. 



i?7r = 7r~/7r+, then its difference from the A isobar 
model prediction i?isob can be used to measure effects 
of the isospin fractionation on pion production. To the 
first order of approximation, the i?^ — i?isob is due to 
the isovector nuclear mean-field. Effects of the latter are 
expected to depend on both the space-time volume and 
the density gradients reached in the reaction. Shown 
in the upper window of Fig. [3] is the mass Agys depen- 
dence of the i?7r — i?isob at beam energies of 0.25, 0.4 
and 0.6 GeV/nucleon , respectively. The results are ob- 
tained using the i?sym(p) with x ~ 1. The i?,r — ^isob 
exhibits a nearly linear dependence on the system mass 
or volume at a given beam energy. At lower beam en- 
ergies, while the maximum density reached is lower the 
reaction time is longer. The net isospin fractionation 
effect of the isovector potential on the 7r~/7r+ ratio is 



thus larger. On the other hand, at higher energies, the 
Rtt — i?isob increases with Agys much faster due to both 
the higher density gradients reached and the larger re- 
action volume available in these reactions. To quantify 
the sensitivity of the -k^ j-n^ ratio to the variation of the 
Esymip) and present its systematics compactly, we de- 
fine the sensitivity as the double ratio of the 7r~/7r+ ob- 
tained with the Esymip) of x = 1 over that with a; = 0. 
The sensitivity is shown in the lower panel of Figl3] as 
a function of the i?^ — i?isob- It is interesting to see 
that at a given value of the i?^ — i?isob, the sensitiv- 
ity is independent of the reactions systems considered as 
they all have the same initial neutron/proton ratio N/Z. 
Moreover, the sensitivity increases with the R^^ — i?isob 
as one expects since effects of the isovector potential is 
approximately linearly proportional to the isospin asym- 
metry of the medium. Overall, all indications are that 
the 7r~/7r+ ratio is probing the strength of the isovec- 
tor potential and depends on the space-time volume of 
the reaction. More quantitatively, for the reactions con- 
sidered, the sensitivity increases from about several per- 
ccnts to approximately twenty percents at the maximum 
isospin fractionation. It suggests that the tt^/tt"*" ratio 
in heavy systems at relatively low beam energies is pre- 
ferred for probing the HD behavior of the i?sym(/o) • 

Having established that the tt^/tt^ ratio depends on 
the isovector mean-field through the isospin fractiona- 
tion mechanism, we now examine more closely how the 
latter works and its dependence on the i?sym(p) • To 
quantify the isospin asymmetry of an excited system dur- 
ing heavy-ion reactions, we define nuko = n + fA*^ + 
iA+ + A-,piikc = P + |A+ + iA" + A++. As a typ- 
ical example, the average ratio < N/Z >= nuko/piike 
at the instant of 8, 10 and 12 fm/c when the maxi- 
mum compression is reached in the reaction of Ca-f Ca, 
Sn+Sn and Au-|-Au at 0.4 GeV/nucleon , respectively, 
is shown as a function of the local density in Fig. [H It 
is seen that the maximum density reached is the highest 
for the Au+Au reaction as one expects. To understand 
the results more easily, we recall that for two regions in 
thermal-chemical equilibrium at densities pi and p2 , and 
isospin asymmetries 6i and 62, one has approximately 
[26t Si ■ Esymi.Pi) = 62 ■ EsymiP2)- While the reaction 
system may not be in complete thermal-chemical equi- 
librium, the above relation and the i?sym(p) shown in 
Fig. [T] help us understand easily the observed correlation 
between the < N/Z > and the local density. Compar- 
ing the < N/Z > obtained with a; = 1 and x = 0, it 
is seen that the supra-saturation density region is more 
neutron-rich with a; = 1 as one expects based on the 
Esymip) shown in figure [TJ On the contrary, the subsat- 
uration density region is more neutron-rich with x = 0. 
The effect from both the sub- and supra-saturation den- 
sity behavior of i?sym (p) will compete and contribute to 
the final observable. As shown in figure [2l the fact that 
the TT"/?!-"*" ratio is higher with x = 1 indicates that 
it indeed reflects the < N/Z > of the supra-saturation 
density region, otherwise the final 7r~/7r+ ratio would 
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FIG. 3; (Color online) (a) the system size and beam energy 
dependence of the degree of isospin fractionation denoted by 
Rtt — -Risob (see text), (b) correlation between the sensitivity 
of the n~ /tv'^ ratio to the Eayui{p) and the degree of isospin 
fractionation. The dashed line is for guiding the eyes. 
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be higher with x = if it was mainly determined by the 
-E'sym(p) at subsaturation densities. On the other hand, 
it is seen in figure [2] that the 7r^/7r+ ratios calculated 
with a; = are also above the predictions of the isobar 
model at low beam energies although the < N/Z > at 
supra-saturation densities is slightly below the N/Z value 
of the reaction system as shown in figure |4l This might 
indicate significant changes in the local N/Z ratios dur- 
ing the time evolution of the reaction since the results 
in figure |4] represent only the density dependence of the 
< N/Z > at the respective instants of the maximum 
compression, not for the entire duration of the reaction. 

Since pions always subject to strong final state inter- 
actions the 7r~/7r+ ratio may also be affected by the 
Esym{p) at sub-saturation densities. To quantitatively 
evaluate effects of the low density Egymip) , let's define 
xl{xh) as the x parameter at sub (supra)-saturation 
densities, effects of the i?sym(p) at different densities can 
then be revealed by comparing calculations using differ- 
ent xl{xh) parameters. As an illustration, shown in fig- 
ureOis the 'k~ jix^ ratio at the chemical freeze-out in the 
Au+Au collisions at 0.4 GcV/nuclcon with different 



FIG. 4: (Color online) The evolution of the {NjZ) as a func- 
tion of the local density at the most compressed stage in 
"^Ca+'^^Ca, ^^■'Sn +^2*Sn and ^^'^Au-H^^'^Au collisions at 0.4 
GeV/nucleon with x = 1 (solid lines) and (dashed lines), 
respectively. 



and xh combinations. Comparing the middle two pan- 
els, it is seen that the saturated value of the (7r^/7r+) 
ratio increases by about 7% when the x parameter is 
changed from —2 (stiff) to 1 (soft) at all densities. In- 
terestingly, if we fix the xh cither at —2 or 1 but vary 
the xl between 1 and —2, the 7r^/7r+ ratio changes by 
less than 2%. On the contrary, if the cc^ is fixed, the 
variation of the xh leads to about a 7% change similar 
to the effect observed when the x parameter is varied at 
all densities also between —2 and 1. Thus, again, the 
Tr~ /n'^ ratio really reflects mainly the Esymip) at supra- 
saturation densities with only very minor influences from 
the Egymip) at sub-saturation densities. 

To further investigate the reliability of the tt" / 7r+ ra- 
tio in probing the HD Esymif) , we now turn to the 
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FIG. 5: (Color online) The 7r~/7r"'" ratio with varying 
Esymip) at sub- (xl) and supra-saturation (xh) densities for 
the head-on collisions of Au+Au at 400 MeV/nucleon. 
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FIG. 6: (Color online) The distribution of As as a function 
of density (a) and time (b), (c) the mean local density of As 
as a function of time. The curve in (c) presents the baryon 
density in the central region as a function of time. 



density profile and time evolution of the produced A res- 
onances through both N + N ^ N + A and tt + iV ^ A 
processes. Shown in Figure [H]Ja) is the total multiplic- 
ity Na of primordial A resonances of all charge states 
as a function of the local baryon density at which the 
A resonances are produced during the entire reaction 
of Au-I-Au at 400 MeV/nucleon with x = (sohd 
square) and x = 1 (open triangle). It is seen that the 
A^A has a major peak around 2po and a minor bump 



around O.lpo- The latter has contributions from both 
final state tt + ^ A resonances and the first chance 
N + N iV + A collisions when the surfaces of the 
two colliding nuclei just start overlapping. Interestingly, 
the majority of A resonances are produced in the high 
density region. Moreover, the A^a is appreciably higher 



with the soft £'sym(p) 



1) than that with the harder 



(p) {x ~ 0) in the high density region. This is con- 
sistent with the total pion multiplicity studied in ref. 

We explore next the dynamics of resonance pro- 
duction. As an example, figure [H] (b) shows the Na av- 
eraged over the entire reaction volume as a function of 
time for the Au-|-Au reaction at 400 MeV/nucleon with 
X = 1. Figure [5] (c) depicts the time evolution of the 
mean local baryon density where the A resonances are 
produced. As a reference the average baryon density in 
the most central cell of 1 fm'^ of the colliding system 
is also shown using the solid curve. By comparing the 
results shown in Fig. [5] (b) and (c), one can see clearly 
that the density where the As are produced correlates 
closely to the central baryon density in the colliding sys- 
tem. More quantitatively, the A'^a reaches its maximum 
at about 20 fm/c when the colliding system is still in the 
compressed phase. Indeed, most of the A resonances are 
produced in the supra-saturation density region. Thus, 
the final 7r~/7r+ ratio is more sensitive to the supra- 
saturation density Egymip) instead of the subsaturation 
one. All these demonstrate again that the 7r^/7r+ ra- 
tio carries effectively useful information about the HD 
Esymip) ■ Nevertheless, it is worth noting that at the 
very beginning of the collision (about 5fm/c) when the 
neutron skins of the projectile and the target penetrate 
each other, the pion-like ratio (7r~/7r"'")iike , defined as 
{tt- +A- + ^A°)/{n+ + A++ + ^A+) which naturally be- 
comes the iT~ /iT~^ ratio at the chemical freeze-out, does 
reaches a peak higher than 3 due to the more abun- 
dant nn collisions leading to more tt~ productions [l^], 
but the actual yield of pions there is very small. Some 
of these pions may escape earlier along the transverse 
direction, they may thus affect significantly the high en- 
ergy tail of the pion squeezed-out. However, they have 
only a very small effect on the final pion ratio integrated 
over all space-time. It is also necessary to mention that 
the study presented here uses the free-space NN inelastic 
cross sections and no mean-field effect for pions is taken 
into account. We notice that in- medium effects on the 
NN inelastic collisions [28l . [29| and pions [sO] may affect 
the pion ratio and they deserve further studies. 

In summary, we carried out simulations of the head-on 
collisions of 48Ca-H48cg^^i24gj^ +^^'^811 and i^Uu+^^Uu 
at beam energies from 0.25 to 0.6 GeV/nuclcon using the 
IBUU04 transport model. The tt~ /tt^ ratio increases 
with increasing the system size or decreasing the beam 
energies, indicating the isovector mean-field is at work. 
A comprehensive analysis of the isospin fractionation 
in these reactions having approximately the same neu- 
tron/proton ratio but different masses, the time evolu- 
tion and the spatial distributions of the As demonstrate 
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clearly that the tt"/??"*" is an effective probe of the HD 
behavior of the nuclear symmetry energy. 
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